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ABSTRACT 

An  approximate  solution  for  the  temperature  distribution 
in  a  thin  circular  disk  having  variable  thermal  properties, 
variable  thickness,  and  a  variable  surface  heat  transfer  co- 
efficient is  presented.   The  method  used  was  to  cast  the 
governing  differential  equation  into  the  equivalent  varia- 
tional problem  and  solve  using  the  direct  method  of  Rayleigh 
Ritz .   Agreement  between  the  approximate  solution  and  an  exact 
solution  for  the  case  of  a  constant  thickness  disk  with  con- 
stant properties  was  excellent.   A  simple  computer  program 
which  may  be  used  to  evaluate  the  temperature  profiles  is 
provided. 
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TEMPERATURE  DISTRIBUTIONS  IN  THIN  CIRCULAR  DISKS  WITH  VARIABLE 
THICKNESS,  THERMAL  CONDUCTIVITY,  AND  HEAT  TRANSFER  COEFFICIENT 
(A  DIRECT  VARIATIONAL  APPROACH  USING  RAYLEIGH  RITZ  APPROXIMATION) 

Applications  requiring  knowledge  of  temperature  distribu- 
tions in  disks  are  many;  probably  the  most  common  being  turbine 
rotors.   Exact  solution  for  the  temperature  distributions  in 
disks  with  variable  properties  is  at  present  impossible  except 
for  a  few  special  cases,  however  approximate  techniques  may  be 
used  which  allow  solution  of  these  problems.   The  very  powerful 
variational  calculus  approach  used  in  this  report  is  one  of 
these  techniques.   The  excellent  agreement  between  the  approxi- 
mate solution  and  an  exact  solution  for  a  particular  case  en- 
courages confidence  in  the  analysis. 

The  primary  motivation  for  this  analysis  was  the  thermal 
stress  problem  in  turbine  rotors  which  requires  knowledge  of 
the  temperature  distribution  (for  the  uncoupled  problem).  Large 
temperature  gradients  may  be  present  which  result  in  large  prop- 
erty variation.   Flexibility  in  the  present  solution  allows  an 
iterative  approach  to  the  problem  of  determining  a  temperature 
distribution  for  this  complex  case.   The  temperature  distribu- 
tion may  be  determined  for  an  arbitrary  distribution  of  prop- 
ertieSj then  the  new  temperature  distribution  used  to  obtain  a 
new  set  of  property  variation  and  so  forth  to  the  desired  de- 
gree of  convergence.   It  is  expected  that  convergence  would  be 
quite  rapid.   A  second  intriguing  application  for  those  who 
enjoy  treading  upon  thin  ice  would  be  to  experimentally  deter- 
mine the  temperature  distribution  in  the  disk  (rotating  or  ndt) 


then  deduce  the  required  heat  transfer  coefficient  distribution 
to  produce  that  temperature  distribution.   This  approach  would 
provide  simple  instrumentation  for  rotating  disks  with  complex 
flow  conditions  existing  on  the  faces  and  still  give  reasonable 
guantitative  heat  transfer  information. 
Formulation  of  the  Problem. 

Consider  an  axisymmetric  thin  circular  disk  of  outside 
radius  R  and  arbitrary  thickness  distribution  b(r).   The  only 
requirement  which  we  place  is  that  the  disk  be  thin  enough   to 
ignore  temperature  variations  in  the  axial  direction.   This 
allows  us  to  treat  the  problem  as  a  normal  heat  conduction 
problem  with  a  uniformly  distributed  heat  sink.   We  may  write 
the  governing  equation: 

Avzr  -Qs^o 

The  heat  sink  term  is  expressed  as  a  function  of  the  heat  trans- 
fer coefficient  on  the  face  of  the  disk. 

2h  IT-Tcd) 
Qs   =   E 

If  the  dimensionless  variables  g , ^  are  introduced  where  T   is 

the  film,  (or  bulk)  temperature  for  the  heat  transfer  coeffi- 

cientj 

e=T-Too  v-r 

then: 

VZ9   -1*9=0 

In  cylindrical  polar  coordinates,  for  axisymmetrical  cases 
and  thin  disks: 
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Define  :    *\  /  *    Z  h  R 

A^)=  IT 

The  boundary  conditions  of  interest  are: 

j^  =  0  at  Y^=-0    J    e=90   at  ix=  I 

Another  set  of  boundary  conditions  would  perhaps  specify  the 
heat  flux  at*w=  1  however  for  two  reasons  this  was  not  done 
in  the  present  study,  first  the  heat  transfer  coefficient  is 
high  for  flow  through  a  turbine  as  a  result  of  high  velocities 
hence  the  wall  temperature  will  closely  approximate  the  bulk 
fluid  temperature;  secondly  it  is  difficult  to  determine  the 
actual  heat  transfer  coefficient  for  the  flow  through  the  turbine 
because  of  the  complex  geometry.   It  should  be  noted  that  the 
problem  of  specified  heat  flux  at  rotor  tip  can  be  handled  in 
a  way  similar  to  the  present  problem. 

The  present  analysis  utilized  a  general  second  degree 
polynomial  to  approximate  the  distribution  of  A(/|)  •   This  was 
written:      A  -, 

^Cvp  =  Arf  +?>\   +  C 

This  distribution  was  chosen  because  of  the  general  ease 
with  which  a  parabola  may  be  fitted  to  a  reasonably  smooth  set 
of  data;  certainly  property  data  will  show  few  discontinuities. 
Notice  that  both  constant  values  and  linear  distributions  are 
subsets  of  this  polynomial.   More  complex  polynomials  may  be 
used  at  a  correspondent  increase  in  algebraic  labor. 

The  final  differential  equation  to  be  solved  takes  the  form: 

-  d  (  „d& 


S[to$  +  (Aif+*n*OnB«o 
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I  attempted  to  cast  this  equation  into  some  form  which  would 
allow  direct  analytic  solution  through  the  various  theorems 
available  for  Bessel  type  differential  equations,  none  of  the 
theorems  used  were  applicable  to  this  form.   This  is  the  dif- 
ferential equation  which  is  solved  by  the  approximate  methods 
of  the  variational  calculus. 
Formulation  of  the  Variational  Problem. 

Allow  us  to  consider  only  the  problem  in  which  the  approx- 
imate solution  will  satisfy  the  boundary  conditions,  we  may 
then  write  the  functional  for  a  homogeneous  differential  equation 
of  the  second  order  as 


1  — i 


For  treatment  of  more  general  cases  which  allow  more  freedom 
in  the  choice  of  the  approximating  function,  the  reader  is  re- 
ferred to  Collatz  "Numerical  Solutions  to  Differential  Equations" 

The  Rayleigh  Ritz  method  may  be  used  to  write  the  solution 
in  the  general  form:  ^ 

Several  specific  polynomials  were  used  in  the  selection  of  the 
most  applicable  one.   These  polynomials  are  listed  in  Table  I 
and  the  resultant  temperature  distributions  are  shown  by  com- 
parison in  Figure  2. 

The  general  approach  to  solution  was  to  substitute  the 
chosen  approximation  into  the  functional  then  carry  out  the 
indicated  integration  to  obtain  a  function  in  terms  of  the 
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Figure  2:   COMPARISON  OF  TEMPERATURE  DISTRIBUTIONS 
USING  VARIOUS  APPROXIMATING  POLYNOMIALS. 
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constants  & .     .   These  constants  were  then  treated  as  variables 

1 

of  the  functional  and  the  extremum  obtained  in  the  usual  way. 

Exact  Solution  for  a  Constant  Thickness  Disk  with  Constant 
Properties . 

A  check  of  the  general  validity  of  the  proposed  approxi- 
mate method  was  made  by  comparison  to  the  exact  solution  of  a 
constant  thickness  disk  with  constant  properties.  Recall  the 
applicable  differential  equation: 


anz     n  an.      ~kb~  e  *~  ° 

with  boundary  conditions: 

for  the  case  ~*i h    -     =   const  =  \ 
The  simple  solution  results  that: 

9   _   XoUX^ 

I   =  hyperbolic  Bessel  function  of  order  zero 
o      r 

This  exact  solution  was  used  as  a  criterion  for  the  accuracy 
of  the  approximation.   As  will  be  seen  the  approximate  technique 
reproduced  this  exact  solution  with  no  distinguishable  error. 
Discussion 

The  various  polynomials  used  are  shown  below     table  I. 

i.    e  •  ©cL  i  +  °i  (.i-n*)l 

4.  e  -e0D  +  <Ui-nl)4f30-n4Y] 

5.    B    =  ©o[l  +  °U>-vf') +-/?(' -^3  v-  pCi-^^J 
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Notice  that  (1)  is  just  the  first  approximation  to  all  the 
other  functions.   By  far  the  best  agreement  to  the  exact  sol- 
ution was  obtained  by  number  five.   It  is  interesting  to  note 
that  numbers  1  and  2  do  not  allow  the  function  much  freedom 
at  the  point    =  0,  while  3  and  4  allow  c*  and/3  to  be  arbitrary 
at  "7/  =  0,  and  5  allows  (A  ,/2>   ,   and  u     to  be  arbitrary  at  77  =  0. 
This  seems  to  indicate  that  the  greater  the  freedom  given  the 
coefficients  the  more  rapidly  will  the  solution  approach  the 
exact.   Polynomial  number  3  allows  for  an  inversion  in  the 
temperature  profile,  a  condition  which  is  not  physically  rea- 
sonable for  the  particular  problem  being  solved.   See  Figure  2 
for  comparison  of  the  various  approximation  functions. 

A  detailed  derivation  for  the  case  number  5  is  presented 
in  the  appendix  for  the  case  of  a  parabolic  distribution  of 
the  variable  A-Cy]")  .   This  was  treated  in  the  form: 

A,  B.  and  C  would  be  specified  by  the  empirical  fit  of  the  ex- 
perimental data  for  a  second  degree  polynomial  for  the  particu- 
lar problem  being  studied.   It  would  be  expected  that  the  dis- 
tribution would  be  nearly  linear  for  most  practical  cases  and 
as  a  result  the  squared  term  coefficient  would  be  small.   An 
actual  turbine  might  be  expected  to  decrease  in  width  as  the 
radius  increased,  the  thermal  conductivity  may  either  decrease 
or  increase  with  radius  (temperature)  depending  upon  the  ma- 
terial, and  the  heat  transfer  coefficient  would  be  expected  to 
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increase  with  radius  due  to  the  higher  relative  velocities. 
If  the  term  ^  is  examined  it  may  be  seen  how  these  properties 
influence  the  parameter. 

A  -     kb 

y^  tends  to  increase  with  radius  and  increased  heat  transfer 
coeff icienfj and  increases  with  decreased  thermal  conductivity 
and  thickness.  A  logical  calculation  procedure  is  to  assume 
the  properties  constant  for  the  first  approximation  then  use 
the  deduced  temperature  distribution  for  further  approximations 
The  results  of  several  computations  for  ranges  of  values  of  A, 
B,  and  C  are  presented  in  Table  II.  This  table  presents  the 
coefficients  for  the  expression  for 

and  also  presents  the  coefficients^,  /3  ,   and  u    f or  the  temper- 
ature distribution: 

%o  =  /   +  ML  i-yj2)  4-/3  (<-  >?4)  4-  p  O-l') 

The  computer  program  also  calculates  the  temperature  distribu- 
tion but  this  vast  printout  would  fill  a  book  and  the  coeffi- 
cients presented  here  allow  for  the  same  information. 
Example  of  Table  useage. 

The  use  of  the  Table  is  completely  straightforward.   Sup- 
pose we  fit  a  parabola  to  \    which  as  the  form: 

\=  .01^  +  ,lr|   +  I 
then  we  enter  the  table  at  the  row  where  A  =  . 01 ,  B  =  .1,  and 
C  =  1.  and  read  from  the  same  row  d   =  .19829, fi   -    .020819,  and 
q     -    .000625.   The  temperature  distribution  is  then: 

33 
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RESULTS  OF  ANALYSIS  FOR  TEMP 
THIN  CIRCULAR  DISKS 
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-.11279E+0Q 
-.19 8 2  9  E  +  0  0 

-.3:    +  0  0 

- .  ■•'  . 

-  )  E  +  0  0 
-.55214E-01 

-.12993E+00 
- . 20563E+0Q 
-.31316E+00 

-  .  4  4  3  6  9  E  +  0  0 

"  +  00 
»,17616E+00 
-.2026SE+00 

-  .23G86E+00 
-.2  1+0  0 
-,33  .  >8E+0  0 
-.34O30E+00 
-.24241E-01 

-  .   6  3  7  E  +  0  0 

6208E  +  0  0 

-  .45441E  +  00 

-.24275E-C1 
.  1 0 1 0  E 


-  0  2 
.12225E-Q1 

.20397E-01 

530G+Q0 
.25537* 
.84348E-01 
.826-   -01 

.  1  0  1  c  i 

,i: 

.25562E+00 

.33542E+Q0 
.319335+00 

.30183E+QG 
.27:,3c;J->00 
.  JE+OO 
.1296. 
.76138E- 
.  4  0  0  9  0  E  -  . 
.12!  . 

.14185E+00 

.2550:  : 

-03 
.4096' 
.12894E- 
.   2 '   -01 
.14188E+QQ 
.25501E+0Q 

-02 


-02 
1E-01 

1E-01- 
QE  +  QQ 

•0  0 

-     . 

-01 
9E-01 

-01 


387 

.  i; 

.i 

.2. 

.1: 

.1270 
.2081 
.  4  6 
.14 
.2552 
.   ' 

:,0  77E-01 


.1012 

37 

554 

.  3353 
192 
.'17 

.2732 

.2146 
293 
,93 

.  .710 
.4335 
.1426 

.6:. 

.8993 


•01 

OE  +  00 
6E  +  00 
1E  +  0Q 

■.  0  0 

•0  0 

■0  0 

•00 

0  0 

■  -  0  1 

-GO 
1E  + 


.19486E- 

.1477 

.69     .-03 

-.2531 
-■.11298E- 
.1--.     -01 
♦11358E-01 

-.23153E- 

-.36  .   -01 

-.1377 

■  .. 

-.10967E- 

-.12819E+0G 

-.26147E 
-.42296E+GQ 
-.91932E- . 5 

-  .  6  9  4  8  2  E  -  : 
-.43085E-Q3 
-.25137E- 
-.2- 

-.1105  .  : 
.  i 

-.73 

-.41989E-03 
- .23107E-02 
- .24226E-01 

-  .11062E+0  0 
.  1 8  9  8  8  E  - 
.6715 

-  ,  3  2 1 6 1 E  ■ 
-.24042E-02 

-  .2  . 
-.110863+00 

.19345E- 
. 14 35 2  3  - 
2.533E- 
-.22463E-02 
-.2554 
-.113263+00 
.•I5i  -01 

.1127:     - 
• 
-.2 

6  7  E  -  0 1 

-  .  1 3  8  0  7  E 
-.9451 
-.10995E+00 

i  +  0  0 
.  +  00 

••  .  . .  '  +  oo 

-  .  1  0 

-.3 

■ 

-  .  1  +00 
-.81 
-.45631E- 
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.100 

.001 

1 

.000 

-.19638E+00 

-.17182E-01 

-.11153E-02 

.100 

.001 

2 

.000 

-.317696*00 

-.43420E-01 

-.372266-02 

•  100 

.001 

5 

.000 

-.462  0  66+0  0 

-.142736+00 

-.264376-01 

.100 

.001 

10 

.000 

-.45439E+00 

-.253616+00 

-.11361E+0C 

•  100 

.010 

.100 

-.24583E-01 

-.71150E-02 

.85506E-04 

.100 

.010 

.500 

-.11029E+0Q 

-.977726-02 

-.32025E-03 

.100 

.010 

1 

.000 

-.19647E+00 

-.17899E-01 

-.10213E-02 

.100 

.010 

2 

.000 

-.31764E+00 

-.440216-01 

-.36993E-02 

.100 

.010 

5 

.000 

-.46187E+00 

-.14304E+00 

-.265575-01 

.100 

.010 

10 

.000 

-.45423E+00 

-.253636+00 

-.113866+00 

.100 

.100 

.100 

-.27626E-01 

-.153196-01 

.17755E-02 

.100 

•  100 

.500 

-.11221E+00 

-.17486E-01 

.99955E-03 

.100 

.100 

1 

.000 

-.19732E+00 

-.24998E-01 

-.116005-03 

.100 

.100 

2 

.000 

-.31711E+00 

-.49963E-01 

-.34931E-02 

.100 

•  100 

5 

.000 

-  .4600  IE- 00 

-.146186+00 

-.27770E-01 

.100 

.100 

10 

.000 

-.45260E+00 

-.253806+00 

-.11629E+00 

.100 

1 

.000 

.100 

-.54991E-01 

-.88956E-01 

•14416E-01 

.100 

1 

.000 

.500 

-.12935E+00 

•  -.868426-01 

.10402E-01 

.100 

1 

.000 

1 

.000 

-.20471E+00 

-.88997E-01 

.56250E-02 

.100 

1 

.000 

2 

.000 

-.31187E+00 

-.103636+00 

-.40053E-02 

.100 

1 

.000 

5 

.000 

-.44227E+00 

-.174076+00 

-.412416-01 

.100 

1 

.000 

10 

.000 

-.43721E+ 00 

-.253835+00 

-  .14112E  +  00 

.100 

10 

.000 

.100 

-.17606E+00 

-.334186+00 

-.97287E-01 

.100 

10 

.000 

.500 

-.20251E+00 

-.317996+00 

-  .  1 1 2  7  6  E  +  0  0 

.100 

10 

■  000 

1 

.000 

-.23062E+00 

-.300376+00 

-.13132E+00 

.100 

10 

.000 

2 

.000 

-.27355E+00 

-.271696+00 

- .16656E+00 

.100 

10 

.000 

5 

.000 

- .33596E+00 

-.21243E+0Q 

-.264905+00 

.100 

10 

■  000 

10 

.000 

-.34827E+00 

-.126405+00 

-.42659E+00 

1 

.000 

.000 

' 

.100 

-.23378E-01 

-.55900E-01 

-.29156E-02 

1 

.000 

.000 

.500 

-.10490E+00 

-.537645-01 

-.58966E-02 

1 

.000 

.000 

1 

.000 

- .18731E+00 

-.56292E-01 

-.945236-02 

1 

.000 

.000 

2 

.000 

-.30418E+00 

-.724796-01 

-.16872E-01 

1 

.000 

.000 

5 

.000 

-.446776+00 

-.149476+00 

-.49011E-01 

1 

.000 

.000 

10 

.000 

- .44492E+00 

-.238846+00 

-.143686+00 

1 

000 

001 

.100 

-.23411E-01 

-.559796-01 

-.29020E-02 

1 

000 

.001 

.500 

-.10492E+00 

-.538386-01 

-.588655-02 

1 

000 

.001 

1 

.000 

-.18732E+00 

-.56362E-01 

-.944626-02 

1. 

000 

001 

2 

000 

-.30417E+00 

-.725376-01 

-.16873E-01 

1 

000 

.001 

5 

000 

-.44676E+00 

-.149506+00 

-.49026E-01 

1 

000 

001 

10 

.000 

-.444906+00 

-.238846+00 

-.143715+00 

1 

000 

010 

100 

-.23708E-01 

-.566936-01 

-.27800E-02 

1 

000 

010 

500 

-.105116+00 

-.545126-01 

-.57961E-02 

1 

000 

010 

1 

000 

-.187416+00 

-.56985E-01 

-.93916E-02 

1. 

000 

010 

2 

000 

- .304136+00 

-.73061E-01 

-.16880E-01 

1. 

000 

010 

5 

000 

- .446586+00 

-.149776+00 

- .49161E-01 

1 

000 

010 

10 

000 

-.4  4  4-75E  +  0  0 

-.238846+00 

-.14395E+00 

1. 

000 

100 

100 

-.26650E-01 

-.637535-01 

-.15997E-02 

1. 

000 

100 

500 

-.10701E+00 

-.61179E-01 

-.49280E-02 

1 

000 

100 

1 

000 

-.18831E+00 

-.63151E-01 

-.88777E-02 

1. 

000 

100 

2. 

000 

-.30374E+00 

-.78248E-01 

-.16973E-C1 

1. 

000 

100 

5. 

000 

-.44488E+00 

-.152476+00 

-.50522E-01 

1. 

000 

100 

10. 

000 

-.44325E+00 

-  .238816  +  00 

- . 14644E+00 

1. 

000 

1. 

000 

100 

-.53228E-01 

-.126996+00 

.64575E-02 

i. 

000 

1. 

000 

500 

-.124096+00 

-.121016+00 

.391836-03 

l. 

000 

1. 

000 

1. 

000 

-.19628E+00 

-.118596+00 

-  .66919E-02 

l. 

000 

1. 

000 

2. 

000 

-.29975E+00 

-.12492E+00 

-.20231E-01 

i. 

000 

1. 

000 

5. 

000 

- .42874E+00 

-.17624E+00 

-.653776-01 

1. 

000 

1. 

000 

10. 

000 

-.42901E+00 

-.23686E+00 

-.17175E+00 

1. 

000 

10. 

000 

, 

100 

-.17543E+00 

-.321846+00 

-.12527E+00 

i- 

000 

10. 

000 

, 

500 

-  .20114E+00 

-.30485E+00 

-  .14106E  +  00 

i. 

000 

10. 

000 

1. 

000 

- .22855E+00 

-.28619E+00 

- .160025+00 

l. 

000 

10. 

000 

2. 

000 

-.270666+00 

-.255396+00 

- .19599E+00 

_ 


1.000 

10 

.000 

5.000 

-.33317E+00 

-.19020E-00 

-■.29612E  +  00 

1.  000 

10 

.000 

10.000 

-.34820E+00 

-.96837E-Q1 

-  .45947E  +  00 

10.  000 

.000 

.100 

-  .  47405E-01 

-.23511E+00 

-.15296E-00 

10  .000 

.000 

.500 

-.98951E-01 

-.21034E+00 

- . 16865E+ 00 

10  .000 

.000 

1.Q00 

-.15326E+00 

-  .  1 8  4  8 1 E  +  0  0 

-.18729E+00 

10  .  000 

.000 

2.000 

-.23557E+00 

-.14715E+00 

-.22079E+00 

10.  000 

.000 

5.000 

-.35688E+00 

-.89189E-01 

-.30881E+00 

10  .000 

.000 

10.000 

-.39249E+00 

-.31690E-01 

-.45201E+00 

10.000 

.001 

.100 

-.47430E-01 

-.23511E+00 

-.15298E+00 

10.000 

.001 

.500 

-.98970E-01 

-.21Q35E+00 

- .16687E+00 

10.000 

.001 

1.000 

-.15328E+00 

-.18482E+00 

-.18731E+0Q 

10.000 

.001 

2.000 

-.23557E+00 

-.14716E+00 

-.22081E+00 

10.000 

.001 

5.000 

-.35688E+00 

-  .89186E-01 

-.30883E+00 

10  .000 

.001 

10.000 

-.39249E+00 

-.31673E-01 

-.45204E+00 

10.000 

.010 

.100 

-.47654E-01 

-.23518E+00 

-.15315E+00 

10.000 

.010 

.500 

-.99140E-01 

-.21042E+00 

-.16904E+00 

10  .  000 

.010 

1.000 

-.15339E+00 

-.18488E+00 

-.18748E+00 

10.000 

.010 

2.000 

-.23561E+00  • 

-.14721E+0Q 

'-.22100E+00 

10.000 

.010 

5.000 

-.35682E+00 

-.891586-01 

-.30907E+00 

10.000 

.010 

10.000 

-.39243E+00 

-.31519E-01 

-.45232E+0Q 

10.000 

-100 

.100 

-.49878E-01 

-.23576E+0Q 

-.15479E+00 

10.000 

.100 

.500 

-.10083E+00 

-.21103E+00 

-.17073E+00 

10.000 

.100 

1.000 

- .15455E+00 

-  .18548E  +  00 

-.18925E-00 

10.  000 

-100 

2.000 

-.23601E+00 

-.14767E+0Q 

-.22291E+00 

10.000 

.100 

5.000 

-.35629'E  +  OO 

-.88864E-01 

-.31139E+QQ 

10.  000 

-100 

10.000 

-.39183E+00 

-.29971E-01 

-.45514E+00 

10.000 

1 

.000 

.100 

-.70669E-01 

-.23929E+00 

-.17246E+00 

10.000 

1 

.000 

.500 

-.11669E+00 

•  -.21487E+00 

- .18881E+0Q 

10.000 

1 

.000 

1.000 

-.16542E+00 

-.18923E+00 

-.20792E+00 

10.000 

1 

.000 

2.000 

-.23980E+00 

-.15020E+00 

-.24288E+00 

10.000 

1 

-000 

5.000 

-.35129E+00 

-.84374E-01 

-.33521E+00 

10.000 

1 

.000 

10.000 

-.38634E+00 

-.13586E-01 

-.48353E*0Q 

10.000 

10 

.000 

.100 

- .19711E+00 

-.12001E+00 

-.42299E+00 

10.000 

10 

000 

.500 

-.21684E+00 

-.98633E-01 

-.44031E+0Q 

10.000 

10 

000 

1.000 

-.23847E+00 

-.74023E-01 

-.46123E+0G 

10.000 

10 

000 

2.000 

-.27338E+00 

-.30394E-01 

-.50112E+00 

10.000 

10 

000 

5.000 

-  .33404E  +  00 

.73913E-01 

-.61142E+00 

10.000 

10- 

000 

10.000 

-.36842E+00 

.21931E+00 

-.78408E+00 

TIME,   0  MINUTES  AND 

51  SECONDS 
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Appendix  A  -  Details  of  Computations 
The  functional  may  be  written: 

Tift  =  j  {hid'f-t  n>(  y  *-&yz0Z  +Cy  f/khj 

Define     r    as   follows: 

The   derivative  becomes: 

.-■ 

\^j    a  /    <V\     -    4  ^      -    kV   V^ 

Squaring 

tV'f  *4oLzv\    +■  IMf  ij4  +  (24^  v  ^p2)  v^  + 

Squaring  the  function. 

'VZ   a  14  lei   tZ8  +-2S1  +  elz  +  f"  +■  V^  +  2  ^  (3   4  2|3&» 

+-2ol5»     -  (2o4  +  Z*L7  +  2<i(3  +2«>t^v\2'  +■ 
(oi?-2(3  -2(*2  -  2=^P  ~2?>v)  r\*    + 
Udl(J    "  2P  -l^-l^1  -ZU^V^  4- 
((^+2e<P)^8     +  2(3l^}10    +  Vz^l 

Define  for  ease  of  algebra. 

Kv  3  |  ^^  ^Z(3  +-2^  +  *c2  4-(32  4-#2  +  Z^p  +2p^  *2d.^ 
K3=   ^-Zf3-Z(3Z-a^p-2(2^ 


K-i  =   ^Z 


Integrals  involving  ^   may  be  evaluated 
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joyVZo/y~jtvj:Az  +  ±*3  +d*+  +&K  +**<> 


The  integral  involving  Cy')     becomes: 

Substituting  these  integrals  the  functional  may  be  written: 

TC'V^o^2    4-1^(3      4-3^  +2(22   +  ^'(Stf     +  3tfZ 

+  (^AHb  +  ^)  K3   +  (&*  +  £b+£-C')K4 
Kit  A  +TTB  +ife£)  Ks"  +CwA  +  TsB  +  ^OKfo 

+-CteA  +  tWb  +tVc~)  K7 

The  functional  is  now  differentiated  with  respect  to  the 
coefficient s  o(, /9,  and  7P   to  obtain  the  extremal.   Define 


for  convenience. 


Li  2  ± A    4  i-  B    4  i 
Lis  5  ikA    ±tt  E  +  roC 


c 


L>2  2  S^A  +^fc+4^ 


>  A  ^  <f  B  +  kc 


Differentiating  with   DC  ; 

SI 


t  ( i  I- 1-1 ,  *  2  tl'2  +  L3)  •*    +  2  ( J  +  L,  -L2  -L3+  U)  9 
+  (3  +2L,  -2L2 -2L4  +2L5)#  +(^1-2^ 
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defining  the  coefficients  of  Qt,  A,    and  $  by  M  j ;  m'  and  m' 
The  derivative  may  be  written:  ' 

g  =  M,'  oL  +  M^  (3    +•  M3  tf    +  M4    =0 

differentiation  with  respect  to  &   yields: 

|j[  =  (f  +  2L,-2L2-2L3  +2L4V  + 

(4  +2L,  -4L3  ^Z^P  +• 

Defining  M" ,  M"   M"  in  a  similar  way: 

J.    2,  3 

finally  differentiation  with  respect  to  *$   yields: 

|^  =  (3+  2L,  -2L2  -  2L4  4-2Ls)^  4- 
s=  +2L,  -2L3  -2Z-4  +2U)Q      + 
(fa  -t-2L,-4L4  ^2L7)!C    -t-(2L/-2Z.^)  =  0 


Tf        " 


Define   M_  ,  M   ,  and  M   in  a  similar  way 

The  system  of  linear  equations  to  be  solved  for  o£ , j3 
and  W   becomes: 


M/"d     *    Al%  (3  4    Atf  P  =  -  At/ 

Once  the  desired  values  for  A,  B,  and  C  are  known  the 
included  computer  program  will  evaluate^,  |9>  and  ^P 
according  to  the  above  method  (see  Appendix  B). 
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APPENDIX  B 
Computer  Program  for  Evaluation  of  Temperature  Profiles 

Program  TEMP  carries  out  the  computations  outlined  in 
Appendix  A.   This  program  is  written  in  Fortran  60  for  use  on 
a  Control  Data  Corporation  1604  digital  computer. 

Coefficients  used  to  describe  the  parameter  "\  in  terms  of 
a  second  degree  polynomial  are  given  the  program  as  input. 

Output  information  consists  of  the  temperature  distributions 
as  functions  of  the  non-dimensional  radius.   The  coefficients 
of  the  temperature  distribution  polynomial  are  also  printed 
out. 

A  fixed  point  variable,  N,  determines  the  number  of  sets 
of  coefficients  to  be  computed.   Symbols  and  formats  for  the 
input  data  are  as  follows: 

Variable  Fortran  Format 

A  A  8F10.5 

B  B  8F10.5 

C  D  8F10.5 

N  110 

A  sample  printout  of  the  program  and  typical  data  are  present- 
ed in  Table  3.   A  sample  output  data  sheet  is  presented  in 
Table  2. 
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Table    3        SAMPLE    PROGRAM   PRINTOUT 


-    9  •  s   I 


-      1     C 


. 


i :  . 

i 


/  / 


- 


z . 


_ 


-    • 


..' 


r 
■-    . 

: 


i 
] 

2  > 

3  j 


- 
=  8  . . 

=-2 

- 

-  4  .  ■ 
3    = 

=  3  •  - 
_     -4» 
=  6.- 


i 


- 
- 


(  : 


2. 


3  i  - 


,      .     D(  I  )  ,  £>•.' 
I  )  /  6  . 

/  1  1  •  H  I  )  /  1  0  • 

1    ;'  • 

i )  /  i  4  • 

-4.       2(3     1-2.        3(1] 

-;  . 
-2.  i 

( I  )-2.*R2(  [  ■     ,*R4( I : 

-  4 .  (       +  2  , 

)  -  2 .        3  (  I  )  -  2  .  ■        (  I  )  +  2  •     <  6  (  I  ] 

1     -  ) 

-2.*       [ I )-2.*R4( I ) +2 

)-2.*R3( I )~2.  )+2. 

- 

:  (  . 


.. 


t  u-i  • ) 

,  -   .    i ; 

,  A  L 
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2t  '•  0RMA1  (31  10.3»3X*3(i£11.5,3X)l 
END 

SUBROUTINE  PARAF1T  (C»8»A) 

this  subroutine  will  fit  a  parabola  to  thrfe  points 

DIMENSION  C<3,3),  5(3),  CM(  3  ,3  )  ,:  A  (3  )  ,  DETCMO)?*  K(3) 
DETC=( C < 1, 1 ) *C( 2,2 }*CC 3  ,3 )+C<]  ,2 ) *C( 2,3 ) *r < 3  ,  ]  )  + 
lCU,3)*C(2,l)*C(.3,2n  -  iCil'9l)*C('2f3)'*Cl3ti)  + 

2C  (  .1  , ,2  )  *C  (  2  ,  1  )  *C  (  3  ,  3  )  +C  (  1 ,  3  )  * C  (  2  ,  2  )  *C  (  3  , 1  )  ) 


1  » 


DO  900 

K  =  K 

DO    vu5  N=l,3 

DO  9  05  J-  1,3 
9  c  5'  CM  (N»J)  "  =  C  (  N  ,  J  ) 

DO  .904  1=1,  3 
904  CM! I ,K)  =  B( I ) 

,rMyf1^irloC?!^1)*CM{1,2,*CM(3,3)  +<-"''i<1.2)*CM(2,3)*CM<3,l)  + 
1  CM ( 1  , 3 ) *CM < 2 , 1 ) *CM ( 3 , 2  )  ) - ( CM ( 1 , l ) *  CM ( 2 , 3 ) *CM (3,2)  + 
2CM ( 1 ,2 ) *CM( 2 , 1 ) *CM (3,3) +CM (1,3) *CM ( 2 , 2 ) -CM ( 3 , 1 ) ) 
900  A(K)  =  DETCM<K)/DETC  ^  U'in 

RETURN 

END 

END 
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